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In 2009, Briane and Milton proved mathematically the existence of three-dimensional isotropic 
metamaterials with a classical Hall coefficient which is negative with respect to that of all of the 
metamaterial constituents. Here, we significantly simplify their blueprint towards an architecture 
composed of only a single constituent material in vacuum/air, which can be seen as a special type of 
porosity. We show that the sign of the Hall voltage is determined by a separation parameter between 
adjacent tori. This qualitative behavior is robust even for only a small number of metamaterial unit 
cells. The combination of simplification and robustness brings experimental verifications of this 
striking sign-inversion into reach. 


It is often stated that measuring the sign of the 
Hall voltage allows for determining whether negatively 
charged electrons or positively charged holes dominate 
the electrical transport in ordinary materials such as 
semiconductors or metals [HHS]. In 2009, Briane and Mil- 
ton [3] proved mathematically that this statement does 
not apply to arbitrary artificial materials: They consid¬ 
ered a three-dimensional simple-cubic arrangement of in¬ 
terlinked nearly touching tori made of a first conductive 
material, a second material in between adjacent tori, and 
a third conductive embedding material (see Fig. 1). They 
showed that the sign of the resulting effective Hall coeffi¬ 
cient Rft can be opposite to that of the Hall coefficients 
of the three constituents Ru, i.e., R • Ru < 0. Further¬ 
more, they proved analytically that the behavior of the 
structure in Fig. 1 is isotropic and that a corresponding 
sign-inversion is not possible in two dimensions 00 . 

Their finding [4] emphasizes that the effective proper¬ 
ties of artificial materials called metamaterials need not 
be intermediate to the properties of their constituents. 
Related examples have served as a catalyzer for the 
emerging field of metamaterials: At some finite angular 
frequency uj > 0 , the effective electric permittivity e e ff (cj) 
of a composite Enam made of, e.g., one isotropic ma¬ 
terial with e{uS) = —1 in vacuum (with e vac (co) = 1 ), can 
assume any value e e ff(cu) G [—00,00]. Likewise, a meta¬ 
material composed of one non-magnetic isotropic con¬ 
stituent with magnetic permeability ju(uf) = 1 in vacuum 
(with p vac (co) = 1), allows for /i e ff(^) £ [—00,00]. By 
combination of these two examples, composites of con¬ 
stituent materials with (real part of the) refractive index 
n(uS) > 0 can lead to n G s(uj) • n(w) < 0 [3j[8j[9]. However, 
for electric current conduction, the possibility of such op¬ 
posite sign has not been known until the work of Briane 
and Milton [4|. Furthermore, their example concerns the 
static case, uo = 0 , whereas the other mentioned examples 
are inherently restricted to finite frequencies uj 7 ^ 0 . 

In this paper, we start by verifying the analytical 
Briane-Milton result using numerical calculations, we 
then drastically simplify their suggestion to a single 
porous material, we show that the sign and magnitude of 
the Hall coefficient can be changed by the geometry of the 
pores, and we demonstrate that the behavior converges 


for an accessible number of metamaterial unit cells. On 
this basis, experimental realizations come into reach. 

In the textbook version of the classical Hall effect illus¬ 
trated in Fig. 2(a), a constant electrical current I x (cor¬ 
responding to an applied voltage U x ) is drawn along the 
x-direction through a cuboid isotropic ordinary material 
with volume L x L y L z . A static magnetic field B z is ap¬ 
plied along the ^-direction. We will use this geometry 
throughout this paper. The Lorentz force leads to a 
potential difference along the //-direction called the Hall 
voltage, which is given by Uu = R^I X B Z /L z , with the 
Hall coefficient Rp = 1/p. The electric-charge density 
p can be positive or negative. Ideally, the Hall voltage 
does not depend on the extent of the sample in the x- 
and //-directions, L x and L y , and scales inversely with 
the sample thickness L z . Thus, one usually considers 
thin films to obtain large Hall voltages, or even a single 
atomic layer in the case of the quantum Hall effect on 
graphene m- Fig. 2(b) illustrates a Hall-bar geometry 
composed of a total of 11x5x1 = 55 metamaterial unit 
cells as to be used in several of the calculations to be 
presented below. 

The magnetic component of the Lorentz force makes 
the configuration effectively anisotropic. A more detailed 
mathematical description thus has to start from the static 
version of the continuity equation for the electric current 
density vector j given by 

v-/= v-Gwo = 0 (i) 


with the scalar electrostatic potential (f = (f(r) and the 
(anisotropic) electric conductivity tensor [3] 
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Here, cjq is the scalar conductivity of the isotropic ma¬ 
terial for B z = 0. For small magnetic fields, i.e., in 
the limit (ctoRh ^) 2 1 , all denominators can be ap¬ 
proximated by 1 and we get the usual proportionality 
Lh oc RhB z . In this tensor form, the mathematics can 
be applied to spatially inhomogeneous materials as well, 
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Figure 2. (a) Scheme of a cuboid isotropic homogeneous ma¬ 
terial exhibiting the classical Hall effect with Hall voltage Un, 
resulting from a constant electric current I x and a perpendic¬ 
ular static magnetic field B z . The shown geometry and coor¬ 
dinate system is used throughout this paper. The small cubes 
with side length a indicate the unit cells of the crystalline ma¬ 
terial with integer numbers of unit cells N x , N yi and N z in 
the x -, y-, and ^-direction. This leads to the sample dimen¬ 
sions L x — N x a , L y m N y a, and L z — N z a. (b) Hall bar 
composed of N x = 11, N y = 5, and N z = 1 metamaterial unit 
cells like shown in Fig. 1(b). The Hall bar contains four metal 
contacts (= equi-potential surfaces) to impose the current I x 
and to pick up the Hall voltage Uu, respectively. 


Figure 1. (a) Crystalline three-dimensional simple-cubic 

metamaterial structure proposed by Briane and Milton [4]. 
(b) One cubic unit cell with lattice constant a. It is com¬ 
posed of three different isotropic constituent materials (the 
tori, the spheres, and the surrounding), all of which have the 
same sign of the Hall coefficient, yet different magnitudes. 
Briane and Milton showed analytically that the resulting ef¬ 
fective metamaterial Hall coefficient is isotropic and can have 
the opposite sign. 


for which a —)> cr(r), with a a(r) and p p(r), and 
hence Rh -T Ra(r) (while B z = const.). 

In this paper, in all calculations for just a single con¬ 
stituent material, we will use constituent material pa¬ 
rameters typical for n-doped silicon (e.g., doped with 
phosphorus) at room temperature [IT], namely <to = 
200 AV _1 m _1 and a doping density (electron density) of 
n e = 10 22 m -3 . This choice leads to a charge density 


of p — — en e , with the electron charge —e = —1.602 x 
10 -19 As, and to a Hall coefficient of Ru = 1/p = 
-624 x 10 6 m 3 A 1 s 1 . These values are not critical 
at all and the results presented below can, in fact, easily 
be scaled to any other parameter combination. 

Our numerical calculations have been performed us¬ 
ing the software package COMSOL Multiphysics (Mumps 
solver). Typically, we discretize the entire Hall bar com¬ 
posed of the metamaterial unit cells and the four metal¬ 
lic contacts (= equi-potential surfaces) by « 2 x 10 6 
tetrahedral elements. All results depicted in this paper 
are in the limit of sufficiently small magnetic fields, for 
which the Hall voltage Uu is proportional to the exter¬ 
nal magnetic field B z . We intentionally avoid, however, 
considering too small magnetic fields for which the off- 
diagonal elements of the above conductivity tensor are 
many orders of magnitude smaller than the diagonal el- 
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Figure 3. (a) Calculated Hall voltage Uu (dots) for the 

Briane-Milton metamaterial shown in Fig. 1 for the geometry 
in Fig. 2(b) versus applied magnetic field B z for N x — 11, 
N y — 5, and N z — 1, i.e., for a total of 55 metama¬ 
terial unit cells. The straight line is a guide to the eye. 
Each constituent alone would show a negative Hall voltage 
in bulk form. Material parameters for the spheres Ru = 
—624 x 10 -6 m 3 A _1 s _1 , ao = 200 AV _1 m _1 , for the tori 
Ru — —3.5 x 10 -11 m 3 A _1 s _1 , ao — 3.5 x 10 7 AV _1 m _1 , and 
for the surrounding Ru = 0m 3 A _1 s _1 and ao = 2 AV _1 m _1 . 
The current is I x — 0.1mA, the geometrical parameters are 
R — 10 /im, r — 1.6 /xm, £ — 1.6/xm (for comparison, £ = 2r 
in Fig. 1), hence a = 4 R — 2(£ + 2 r) — 30.4 /xm. The cylin¬ 
drical pick-up contacts have radius R + r, height 0.2 x R , 
and conductivity 3.5 x 10 7 AV -1 m -1 (aluminum), (b) Cal¬ 
culated Hall voltage (dots) for a bulk material with the same 
constituent-material parameters as the spheres in panel (a). 
Note the different vertical scales. The straight line is a guide 
to the eye. 



ements, leading to numerical difficulties. Furthermore, 
the meshing is the same for each unit cell and symmet¬ 
ric with respect to the ^/-direction. Otherwise, artificial 
transverse “Hall” voltages can occur for zero magnetic 
field. 

We start our discussion by considering the geometry 
and the parameters chosen by Briane and Milton [4] (see 
Fig. 1) with a Hall bar composed of 55 metamaterial unit 



Figure 4. (a) Simplified three-dimensional metamaterial 

(compare Fig. 1) composed of a single constituent ma¬ 
terial only. We take n-doped silicon with conductivity 
ao = 200 AV -1 m -1 and Hall coefficient Ru = —624 x 
10 -6 m 3 A -1 s -1 . (b) One cubic unit cell with lattice con¬ 
stant a. The porous structure can lead to a Hall voltage with 
inverted sign with respect to that of the homogeneous bulk 
material (see Fig. 5). A crucial parameter is the separation d 
between adjacent tori, which can be positive, zero, or negative 
(see Fig. 5). The depicted configuration corresponds to d < 0. 
The torus diameter is 2 R and the torus wire diameter 2r. In 
this paper, we fix R = 10 /xm. The lattice constant a results 
from a — 4 R + 2d. 


cells like shown in Fig. 2(b). In Fig. 1(b), the spheres 
touch the tori at singular points only. This configuration 
is pathological numerically. We have thus rather consid¬ 
ered a small but finite contact area. Corresponding nu¬ 
merical results are depicted in Fig. 3(a). They reproduce 
the analytically predicted sign-inversion of the Hall volt¬ 
age with respect to the bulk material shown in Fig. 3(b). 

Such intricate three-dimensional microstructures con¬ 
taining three different constituents can possibly be re¬ 
alized experimentally, but their microfabrication is quite 
demanding by today’s standards. A central aspect of this 
paper is thus to simplify the microstructures as much as 
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Figure 5. (a) Calculated Hall voltage Uu (dots) for the sim¬ 
plified single-component metamaterial shown in Fig. 4 for the 
geometry in Fig. 2 and for fixed I x — 0.1mA versus ap¬ 
plied magnetic field B z for various parameters. The straight 
lines are guides to the eye. In bulk form, the constituent 
material would lead to a negative Hall voltage, (b) Calcu¬ 
lated Hall voltage (dots) versus separation parameter d (see 
Fig. 4) divided by the torus radius R. The curves are guides 
to the eye. Obviously, the sign of the Hall voltage corre¬ 
lates with the ratio d/R, whereas the wire radius r has a 
much smaller influence. The insets illustrate the geometry 
for three different values of d/R. Parameters are: R = 10 /xm, 
B z = 1 T, N x = 11, N y = 5, N z = 1, contacts as shown in 
Fig. 2(b), and n-doped silicon with ao = 200AV -1 m -1 and 
Rh = -624 x 10 -6 m 3 A -1 s -1 . 


possible, which is possible indeed: We obtain similar re¬ 
sults for using only a single constituent material in vac¬ 
uum/air (illustrated in Fig. 4, n-Si parameters as quoted 
above). For example, for a moderate magnetic field of 
B z = IT and a constant current of I x = 0.1 mA, we cal¬ 
culate a Hall voltage in the range of 10 mV as shown in 
Fig. 5(a). The combination of these parameters appears 
quite reasonable for practical experimental realizations. 

Figure 5(b) depicts results for different separations of 



N x 


Figure 6. Calculated Hall voltage Uu (dots) for fixed ge¬ 
ometry and parameters of the metamaterial unit cell ver¬ 
sus N x for different values of N y as indicated in the legend. 
We fix N z = 1. The shape of the contacts is as shown in 
Fig. 2(b). The other fixed parameters are cr 0 = 200 AV -1 m -1 , 
Rh = -624 x 10 -6 m 3 A _1 s _1 , R = 10/xm, r/R = 0.1, 
d/R — —0.44, B z — IT and I x — 0.1 mA 



Figure 7. Hall potential map, i.e., difference of the calcu¬ 
lated electrostatic potential <f>(r) with [B z = IT) and without 
(B z = 0) magnetic field for the parameters as in Fig. 6 for the 
case of N x = 11, N y = 5, and N z — 1 plotted on a false-color 
scale. 


adjacent tori, d (as defined in Fig. 4), while fixing the 
constituent material and while fixing the radius of the 
tori to R = 10 (jl m. As a result, the metamaterial lattice 
constant a changes according to a = 4 R + 2d (compare 
Fig. 4). Obviously, the sign of the Hall voltage changes 
versus d: For thin tori wires, r/R <C 1 (see Fig. 4), the 
zero crossing of the Hall voltage Ua(d) asymptotically 
occurs at d = 0, whereas it is somewhat shifted from 
zero for larger ratios r/R. 

For bulk Hall media, it is well known meg that the 
actual Hall voltage depends on the geometry. In case of 
cuboid geometry like shown in Fig. 2, the thickness L z 
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Figure 8. Illustration of the sign-inversion of the Hall voltage 
for different tori separations d. From top to bottom: d > 0, 
d ~ 0, and d < 0. The schemes in the left-hand side column 
can be compared with parts of the metamaterial structure 
shown in the right-hand side column. 


is usually chosen to be small compared to the other di¬ 
mensions to maximize the Hall voltage. Furthermore, the 
width L y has to be small compared to the length L x and 
the pick-up contacts have to be small compared to the 
width in order to minimize the influence of the contacts. 
In addition, for metamaterials, unlike for ordinary crys¬ 
tals, the unit cell is no longer very small compared to the 
size of the Hall bar. We have thus investigated the depen¬ 
dence of the Hall voltage on the number of metamaterial 
unit cells for fixed shape and parameters of the unit cell. 
Results are depicted in Fig. 6 . Altogether, these results 
suggest that a Hall bar composed of N x = 11, N y = 5, 
and N z = 1 metamaterial unit cells, under the quoted 
conditions, together with the contacts shown in Fig. 2(b) 
approximates reasonably well a fictitious infinitely ex¬ 
tended Hall-effect continuum. If the pick-up contacts (= 
equi-potential surfaces) are artificially removed, the Hall 
voltage changes by a very few percent with respect to the 
case with pick-up contacts (Fig. 5). 

In Fig. 7, we depict the difference of the electrostatic 
potential </> with ( B z = 1 T) and without (. B z = 0) mag¬ 
netic field for the parameters as in Fig. 6 for the case of 


N x = 11, N y = 5, and N z = 1. Obviously, the edge 
effects decay towards the center of the Hall bar, again in¬ 
dicating that we approach the conceptual ideal of an in¬ 
finitely extended Hall-effect continuum with a finite and 
potentially accessible number of metamaterial unit cells. 
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Figure 9. (a) Hollow version of the single-component meta¬ 
material in Fig. 4. The shown mesh has been used for the 
numerical calculations, (b) Calculated Hall voltage (dots) 
as Fig. 5(b) but for the hollow structure shown in panel (a). 
Parameters are: N x = 11, N y = 5, N z = 1, contacts as 
shown in Fig. 2(b), I x m 0.1mA, B z = 1T, R = 10 /mi, 
r/R = 0.14, and n-doped silicon with ao = 200AV -1 m -1 
and Ru = —624 x 10 -6 m 3 A -1 s -1 . The three different coat¬ 
ing thicknesses t/r are indicated in the legend. 

The dependence Uu(d) shown in Fig. 5(b) is the key 
to an intuitive understanding of the Hall-effect sign- 
inversion in the single-constituent metamaterial. To start 
with, note that the tori parallel to the yz- plane lead to 
negligibly small local Hall voltages along the ^/-direction 
because the vector product of current and magnetic field 
has no ^/-component for these tori. The situation is differ¬ 
ent for the tori parallel to the xy- plane and those parallel 
to the xz- plane. 

Consider the tori parallel to the xy- plane illustrated 
in Fig. 8. At their center along x, the current essentially 
flows in the x-direction. Together with the magnetic field 
parallel to z, this current leads to an ordinary Hall volt¬ 
age along the ^/-direction across the torus wire thickness 
on the left as well as on the right-hand side of the torus. 
For d > 0, the tori parallel to the yz- plane essentially pick 
up this voltage. In contrast, for d ~ 0, these tori connect 
to points of equivalent potential and the external Hall 
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voltage is zero. For d < 0, the yz -torus on the left of 
the xy-toius connects to the right-hand side of the torus 
wire and vice versa. As a result, the yz -tori pick up the 
local Hall voltage and guide it to the outside with an in¬ 
verted sign. This reasoning implies that the separation d x 
along the x-direction is of minor importance—although 
the current locally flows in the negative x-direction at 
these connections. In additional calculations (not de¬ 
picted) we have varied d x ^ d while fixing the separa¬ 
tions d = d y = d z . Indeed, the sign of the effective Hall 
voltage does not change when changing the sign of d x . 

The overall reasoning is analogous for the Hall voltage 
due to the tori parallel to the xz-plane. Again, the tori 
parallel to the yz -plane pick up the local Hall voltage and 
guide it to the outside world. 

In a strictly two-dimensional system parallel to the 
xy- plane, the crucial yz -tori do not exist. Thus, intu¬ 
itively, the Hall-effect sign-inversion is not possible in 
two dimensions—as was previously proven mathemati¬ 
cally [HE]. 

Let us note in passing that our intuitive reasoning sug¬ 
gests that a related sign-inversion is expected for the 
Ettingshausen-Nernst effect. Herein, a temperature gra¬ 
dient along the x-direction (instead of a potential gra¬ 
dient for the Hall effect) leads to a heat current and to 
an electric current in the x-direction. Together with a 
magnetic field along z, in analogy to the Hall effect, this 
current leads to a transverse voltage in the ^/-direction. 
Again, depending on the tori separation d, i.e., on how 
this voltage is picked up internally, we expect a positive 
or negative external Ettingshausen-Nernst voltage—for 
the same constituent material. 

Single-component structures like shown in Fig. 4 could 
be fabricated by three-dimensional direct laser writing 
(3D DLW) p~5] of a polymer scaffold, followed by repli¬ 
cation using a double-inversion procedure for silicon m 
or titania m This technology has originally been devel¬ 
oped for photonic crystal fabrication. Recently, 3D DLW 
has been enhanced by rapid galvo scanners [H> such that 
structures with sub-micron feature sizes and, at the same 
time, overall volumes of some mm 1 2 3 4 5 6 7 8 can be fabricated 
within some hours of writing time. However, double¬ 


inversion procedures are still quite demanding techno¬ 
logically. It is much simpler and more reliable to just 
conformally coat the polymer scaffolds by a semiconduc¬ 
tor material. For example, promising results have been 
obtained by using atomic-layer deposition of ZnO doped 
with A1 (“AZO”) at variable concentrations HMD- The 
resulting geometry, however, is distinct from the one con¬ 
sidered so far. 

We have thus repeated the calculations for hollow 
tori (illustrated in Fig. 9(a)) instead of massive tori (see 
Fig. 4). Hollow tori are a reasonable approximation be¬ 
cause the polymer scaffold can be considered as an iso¬ 
lator compared to the semiconductor and is thus elec¬ 
trically equivalent to vacuum/air. For consistency, we 
continue using parameters for n-doped silicon as con¬ 
stituent material. Corresponding results are exhibited 
in Fig. 9(b). Again, the overall qualitative behavior is 
qualitatively rather similar to what we have discussed 
above, once again highlighting the robustness of the sign- 
inversion effect for the single-constituent metamaterial. 
The moduli of the Hall voltages for r/R = 0.14 are about 
six times larger for the hollow than for the solid tori be¬ 
cause the current is restricted to a thinner layer along 
the ^-direction. As discussed in the introduction, the 
Hall voltage scales inversely with the thickness in the z- 
direction. 

In conclusion, we have shown that the Hall voltage of 
a suitably shaped three-dimensional porous material and 
that of an otherwise identical bulk material can have 
opposite sign. This means that, in effect, a hole con¬ 
ductor appears like an electron conductor or vice versa. 
This finding challenges the common textbook wisdom on 
the connection between the dominant type of charged 
carrier and the sign of the Hall voltage. We have ar¬ 
gued that it should be possible to realize such metama¬ 
terial structures using state-of-the-art three-dimensional 
micro-fabrication technology. On the basis of our intu¬ 
itive discussion of the calculation results, we expect a re¬ 
lated sign-inversion for the Ettingshausen-Nernst effect 
on the same metamaterial structures. 

We thank the Hector Fellow Academy and the Karl¬ 
sruhe School of Optics & Photonics (KSOP) for support. 
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